Evolution of Oscillating Scalar Fields as Dark Energy 



Sourish Dutta and Robert J. Scherrer 
Department of Physics and Astronomy, Vanderbilt University, Nashville, 



TN 37235 



00 

o 

O 

(N 



Oscillating scalar fields, with an oscillation frequency much greater than the expansion rate, have 
been proposed as models for dark energy. We examine these models, with particular emphasis on 
the evolution of the ratio of the oscillation frequency to the expansion rate. We show that this 
ratio always increases with time if the dark energy density declines less rapidly than the background 
energy density. This allows us to classify oscillating dark energy models in terms of the epoch at 
which the oscillation frequency exceeds the expansion rate, which is effectively the time at which 
rapid oscillations begin. There are three basic types of behavior: early oscillation models, in which 
oscillations begin during the matter-dominated era, late oscillation models, in which oscillations 
begin after scalar-field domination, and non-oscillating models. We examine a representative set 
of models (those with power-law potentials) and determine the parameter range giving acceptable 
agreement with the supernova observations. We show that a subset of all three classes of models 
can be consistent with the observational data. 
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I. INTRODUCTION 

It appears that approximately 70% of the energy den- 
sity in the universe is in the form of an exotic, negative- 
pressure component, dubbed dark energy [J 0, 0, 51] ■ (See 
Ref. @ for a recent review). The most likely possibility 
for this dark energy is a cosmological constant. Another 
possibility, dubbed quintessence, is a model in which the 
dark energy arises from a scalar field [1, 0, H, HI E3| • 

Although it is usually assumed that the evolution of 
the scalar field is monotonic, a more exotic possibility 
is that the scalar field oscillates. Such oscillating scalar 
fields were first investigated by Turner (llj . and later 
reexamined in the context of models for dark energy 
lllllllll- (See also Ref. El for a phantom 
version of these models, and Ref. [17| for applications 
to inflation). Here we examine a key aspect of such 
models: the interplay between the oscillation frequency 
and the Hubble expansion rate. The discussions in Refs. 
[ill . El, El, E3, EH assumed that the oscillation frequency 
is large compared to the expansion rate. However, if the 
ratio between these two quantities crosses unity, it is pos- 
sible to generate interesting new behaviors for the dark 
energy. We explore this possibility here. (Note that the 
models considered here, in which the scalar field oscil- 
lates, are distinct from models in which the equation of 
state oscillates 0, El, ES IH SI- In the latter case, 
one is generally interested in oscillation frequencies which 
are always on the order of the Hubble parameter. In the 
models considered here, oscillations with frequencies on 
the order of the Hubble parameter are a transient, albeit 
important, phenomenon). 

In the next section, we review the evolution of oscil- 
lating scalar fields in an expanding background, and we 
derive some interesting general results on the evolution 
of the ratio of scalar field oscillation frequency to the 
expansion rate. In Sec. Ill, we compare numerical simu- 
lations of a representative set of oscillating models (due 
to power-law potentials) with the observations. Our con- 
clusions are summarized in Sec. IV. Our main new result 
is the existence of a class of oscillating models in which 



the scalar field is initially slowly rolling, but oscillates at 
late times; we show that such models can be consistent 
with the observations, but this result is very sensitive to 
the phase of the scalar field oscillation at the present. 



II. EVOLUTION OF OSCILLATING SCALAR 
FIELDS 

A. Previous work 

Consider a minimally-coupled scalar field, </>, in a po- 
tential V(4>). The equation of motion for this field in an 
expanding background is 



3H </> + -77 = °> 



where the Hubble parameter H is given by 



(1) 



(2) 



Here a is the scale factor, px is the total density, and 
we work in units for which 8irG = 1. The pressure and 
density of <f> are given by 



h= - =vVA 



and 



P0 - y + v(4>), 



(3) 



(4) 



respectively, and the equation of state parameter, w^, is 
defined to be 



=p<p/p < f>. 



(5) 



The density of the scalar field then evolves with the scale 
factor as 



P'i> 



(6) 
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Turner was the first to consider the behavior of 
such a field in the limit where it oscillates on a timescale 
much shorter than the Hubble time. He examined power- 
law potentials of the form 



v(<f>) = W, 

and he showed that the value of w,p, 
oscillation period, is 



(7) 

averaged over the 



W4> = (n-2)/(n + 2). 



(8) 



Thus, the cases n — 2 and n — 4 behave like nonrela- 
tivistic matter and radiation, respectively. 

Liddle and Scherrer investigated oscillating scalar 
field models with power-law potentials in the context of 
dark energy, noting an interesting set of behaviors: for 
a universe dominated by a background fluid with equa- 
tion of state parameter wb, the oscillating solution is an 
attractor of the equations of motion only for the case 



n < 



2(3 + w B ) 



(9) 



while potentials with n > 2(3 + w_b)/(1 — wb) produce 
a scalar field trajectory that decreases monotonically to 
zero without oscillating. Liddle and Scherrer provided 
no physical explanation for this curious behavior; we give 
such an explanation below. 

Sahni and Wang [12] and Hsu [HI examined the pos- 
sibility that an oscillating field could give rise to dark 
energy. From equation (JHJ) , it is clear that w < requires 
n < 2, and the production of a dark energy equation of 
state can only be achieved if we take V = k\4>\ n , with 
n -C 1. Although potentials with very small n might ap- 
pear rather unnatural, such a potential can arise from a 
Lagrangian with a <f> 2 potential and a non-standard ki- 
netic term [l^ . 

Masso et al. [l4| performed a general investigation of 
oscillating scalar fields in arbitrary potentials (as opposed 
to power laws). Using the action-angle formalism, they 
derived several very important results, which we will use 
here. Consider a scalar field (f> oscillating in an arbitrary 
potential V ((/)). The minimum and maximum values of 4> 
are <j>\ and (f>2 respectively. Then the energy density for 
this field must be equal to = V{4>\) — V{4>2)- (Note 
that p0 changes with time, but it is assumed to evolve on 
a timescale longer than the inverse oscillation frequency) . 
Masso et al. showed that all of the quantities of interest 
can be expressed in terms of the action variable J, defined 
by 



J = 2 J 2( 94> -V (<!>))■ (10) 

Hi 

Then the oscillation frequency v is 



while the equation of state parameter is 
J 1 



1 



W<j, 



p^ dJ/dp^ 

2 ^[p,~v(4>)]^d4> 



(12) 
(13) 



(Note that this last result was also derived by Turner 
[111]). These results allow for the investigation of some- 
what more natural (non-power-law) potentials that might 
serve as dark energy and several of these are discussed 
inRef. [3. 

Finally, Gu [l5| revisited the power-law scenario, and 
calculated the conditions necessary for the oscillation fre- 
quency to exceed the expansion rate (our equation [16] be- 
low). He used this condition to place constraints on such 
power-law models. 



B. Evolution of oscillation frequency 

Implicit in all of these discussions is the assumption 
that the oscillation frequency is much greater than the 
Hubble expansion rate. We now consider this condition 
in more detail. The relevant quantity is the ratio of oscil- 
lation frequency to Hubble parameter, v/H. If v/H ^> 1, 
then the parameters of interest can be averaged over the 
oscillation period, as in Refs. [|, M, H 51 Q O- 
On the other hand, when v/H <C 1, the field effectively 
ceases to oscillate, and an alternative treatment is re- 
quired. 

Consider the evolution of v/H. From equations (fT0|) 
and pip we can write 



4>i 



P4> ~ V{4>) 



(14) 



This equation, along with the expression for H given by 
equation gives us 



v/H 



2 PT 



3^0 - vm 



(15) 



For the case of the power-law potential, V(<p) = k(f> n , the 
expression for v can be integrated exactly [bj ]. giving 
(see also [l5j ) 



v/H 



V3nT(l/2+l/n) 



2V27rr(l/n) 



(16) 



where 4> m ax = 4>i = 02, since the potential is symmetric 
For other symmetric potentials, we can use equation Q15[ 
to derive an order of magnitude estimate: 



(11) 



(17) 
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This result is clearly consistent with the exact result for 
power-law potentials given by equation (fTB]) . 

Consider first the case where the scalar field itself dom- 
inates the expansion, so that px ~ p^. In this case, we 
have v/H ~ l/4>max- Since p^ always decreases with 
time, 4>max is also a decreasing function of time, so that 
v/H always increases with time. 

The other possibility is that the universe is dom- 
inated by a background fluid (e.g., matter or radia- 
tion) with density ps- In this case, we have v/H ~ 
\J 9$/ PBiX/fymax)- We can again assume that 4> max de- 
creases with time. In order for <j> to serve as a plausible 
dark energy candidate, its density must decrease more 
slowly with the expansion than the density of matter or 
radiation. Hence, we can conclude that 9$/ 9b increases 
with time, and once again we have that v/H increases 
with time. 

Thus, for any model in which <j) serves as the dark en- 
ergy today, v/H is an increasing function of time, regard- 
less of whether the universe is background-dominated or 
^-dominated. 

Before examining the implications of this result, we re- 
consider the behavior discussed in Ref. @ and described 
by equation ([9]). Equations (|5J) and ©, taken together, 
imply that scalar field oscillations are an attractor in a 
universe dominated by a background fluid only for the 
case where 

w <j ,<l/2 + w B /2. (18) 

Further, for a power-law potential given by equation l[7]). 
in the background-dominated case, equation p7[) can be 
expressed as 

ITT -1/2 1/2-1/n ., n x 

V I H ~ P B ■ ( 19 ) 

Applying equation ©, this becomes 

v/H-p-^py'^, (20) 

which can be rewritten in terms of the scale factor (using 
equation © for and the corresponding equation for 
Pb) as 

v/H~a 3 W 2+WB / 2 - w *\ (21) 

A comparison of equations (f2"Tj) and (fT5|) provides a phys- 
ical explanation for the bound noted by Liddle and Scher- 
rer [§]. When this bound (equation [9J is not satisfied, 
equation (|2~lj) indicates that v/H is a decreasing function 
of time. Thus, at late times, the oscillating solution can- 
not be an attractor. This is consistent with the behavior 
noted in Ref. [9(; when equation ([9]) is not satisfied, the 
scalar field decays asymptotically to the minimum of the 
potential without oscillating. These models do not con- 
tradict our claim that in any oscillating scalar field model 
for dark energy, v/H always increases with time. The 
reason is that models that violate the bound in equation 
^ have an energy density that decays more rapidly than 
Pb', they are not plausible models for dark energy. 



Now consider the possible oscillating scalar field mod- 
els for dark energy. The fact that v/H increases with 
time simplifies the classification of such models. If we as- 
sume that the universe evolves from a matter-dominated 
state to a scalar-field dominated state, there are basically 
three possibilities, which depend on the redshift at which 
v becomes greater than H: 

I. The oscillation frequency becomes greater than the 
expansion rate during the matter-dominated era. Since 
v/H increases with time, we always have v/H > 1 during 
the quintessence-dominated era, and we can use time- 
averaged quantities to describe the equation of state of 
the quintessence. We will call this the early-oscillation 
case. This is basically the model described in Refs. [ill, 
El, EE EH EH j so these models have been exhaustively 
explored. 

II. The oscillation frequency is less than H at the 
present day. Again, since v/H always increases with 
time, this implies that v is never larger than H, so these 
models never reach the oscillatory stage. In the language 
of Ref. [23| , these are thawing models in which the field 
has not yet reached the minimum of the potential. Such 
models are characterized by a value of w$ very close to 
— 1 initially, which increases at late times. We will call 
this the "no oscillation" case. Again, such models have 
been exhaustively explored (see, e.g., Ref. [24[ and ref- 
erences therein). 

III. The case intermediate between I and II occurs 
when v becomes greater than H at late times, after the 
quintessence density exceeds the matter density. In this 
case, the scalar field displays two distinct behaviors dur- 
ing the quintessence-dominated era. At first, when the 
oscillation frequency is less than H, the scalar field rolls 
down the potential and w$ increases. Then, when the 
oscillation frequency becomes greater than H, the field 
oscillates rapidly, and t akes its time-averaged value 
as in Refs. [H, El El, M El ■ This third possibility 
has not been previously studied; we will refer to it as the 
"late oscillation" case. 



III. COMPARISON WITH OBSERVATIONS 

Although we believe that our classification scheme in 
the previous section should describe most possible mod- 
els for quintessence from oscillating scalar fields, it is ob- 
viously impossible to compare all such models with the 
observational data. We will therefore limit ourselves to 
a representative set of models: the power-law oscillating 
scalar field models discussed in Refs. [HI, 

ElElEl- As 

we note below, we believe that the general features of the 
results we derive for these models will be more generally 
applicable to other oscillating scalar field models. 

Following Ref. [Hj], we choose the potential to have 
the form: 

V{<t>) = k\4>\ n , (22) 
where k and n are constants. Then the evolution of the 
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field depends on three parameters: k, n, and the initial 
value of fa which we denote fa. However, the parameter 
k can be eliminated from the equation of motion (equa- 
tion [1]) by a suitable rescaling of <fi and the time t. This 
rescaling means that we are no longer using Planck units, 
but our units of time are irrelevant, since we always de- 
fine the "present day" in terms of the value of f^, rather 
then t. (Another way to see this is to note that all of our 
results depend only on the evolution of as a function 
of fi^, and uu[TL,) depends only on V'/V, which is inde- 
pendent of k [IJJ]). Thus, the evolution can be described 
entirely in terms of n and fa. 

For a given choice of values of the two free parame- 
ters, fa, and n, we numerically evolve the system from 
a starting point in the matter dominated regime to the 
point when Q m — 0.3, which we label as the present time 
z = Q. 

The relation between the evolution of <f> and the two in- 
put parameters is illustrated in Fig. [TJ In our numerical 
simulations, we define oscillations to commence when the 
scalar field crosses (f> = 0. This corresponds to w$ increas- 
ing from —1 to +1 and marks the point at which dw^/dz 
first changes sign. (Note that the oscillation frequency 
for is twice as large as the oscillation frequency for 
fa). For early oscillations, this zero-crossing takes place 
during the matter-dominated era, while for late oscilla- 
tions it takes place during the quintessence-dominated 
era. The no oscillation case corresponds to the parame- 
ter region in which zero crossing never occurs up to the 
present. 

For a given n, the parameter fa determines when the 
oscillations will begin. From equation (fTB)) (or [T7|) , it 
is clear that small fa corresponds to a larger value of 
v/H for a given ratio of to px- Thus, oscillations 
commence earlier if fa is small and later if fa is large. 
For fixed n, sufficiently small fa will produce the early 
oscillation case, while sufficiently large fa gives no oscil- 
lations, with late oscillations arising from intermediate 
values of fa. This is apparent in Fig. [TJ Fixing n, we 
see that increasing fa moves the evolution from early os- 
cillation, through late oscillation, and into no oscillation. 
Since the late oscillation case has not been previously ex- 
plored, we present a graph, in Fig. [5J of the evolution of 
w(z) for some of these models. As claimed, initially 
evolves slowly upward from = — 1 , but then oscillates 
rapidly at late times; this transition takes place after the 
quintessence field becomes dominant. 

Next, we compare the theoretical luminosity- vs- 
distance results to the recent Type la Supernovae stan- 
dard candle data (ESSENCE+SNLS+HST from [|). 
The theoretical distance modulus in our model is cal- 
culated in the standard way: 

H{z) = 51og 10 D L {z) + 42.38 - 51og 10 h, (23) 

where h is the Hubble parameter in units of 100 km sec -1 
Mpc -1 , and Dl(z) is the dimensionless luminosity dis- 
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FIG. 1: For the potential V(<j>) = k\<f>\ n , in the parameter 
space determined by n and the initial magnitude of the scalar 
field, fa, this plot gives the regions where the field <j> oscillates 
prior to dominating the expansion (early oscillations), after 
dominating the expansion (late oscillations), or does not os- 
cillate at all prior to the present (no oscillations). 
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FIG. 2: Evolution of the quintessence equation of state pa- 
rameter, w$, as a function of redshift z, in the late oscillation 
regime. Curves correspond to the indicated values of the ini- 
tial value of the field, fa, in the power-law potential V = k\<f>\ n 
with n = 0.16. Diamonds mark the points at which the den- 
sity of dark energy equals the density of matter. The curve 
for fa = 0.15 corresponds to the no oscillation regime and has 
been added for comparison. 



tance, given by 

D L {z) = (1 + z) 



H 



o H(z') 



-dz', 



(24) 



where H(z'), the Hubble parameter at redshift z' , also 
depends on f2 m o and the parameters of our quintessence 
model. We construct a x 2 likelihood plot for fa and n 
ranging from 0.01 to 0.40 (no further structure emerges 
from extending the parameter space). We take fl m o = 0.3 
and marginalize over the nuisance parameter h. Fig. [3] 
shows this plot. 

Some of the allowed regions can be understood in terms 
of our previous discussion of the behavior of the field. In 
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FIG. 3: Likelihood plot for the parameters <f>i and n. The 95% 
confidence region is green (dark shaded) and the 68% confi- 
dence region is green+ yellow (dark shaded + light shaded). 
The diamonds indicate the models displayed in Figs [2] [5] [6] 
and [3 



no oscillations 




FIG. 4: The likelihood contours of Fig. [3] superimposed over 
the demarcation of the type of oscillatory behavior from Fig. 
□ 



the limit of large </>i, the field never oscillates, and it can 
be effectively treated as a cosmological constant, which is 
known to be a good fit to the observational data. For suf- 
ficiently small <j>i, we are in the "early oscillation" regime. 
In this case, takes on a constant value determined by 
equation f8} and independent of 4>i . It is clear from equa- 
tion ([5]) that a value of near —1 requires a value of n 
close to zero. Thus, we see in Fig. [3] that for small </>i, 
the allowed region is confined to small values of n (which 
give u>0 near —1), and the likelihood contours become 
vertical lines (indicating that the late-time value of Wtf, is 
independent of 4>i in this limit). 

To explore the allowed parameter space further, we 
have superimposed the field behaviors outlined in Fig. [1] 
onto the likelihood contours displayed in Fig. [3J The re- 
sults are displayed in Fig. [4j As noted earlier, almost the 
entire "no oscillation" parameter space is allowed, while 
the "early oscillation" parameter space is confined to the 
region for which is close to —1. However, we also 



see that there is an allowed region in the "late oscilla- 
tion" regime; this is the subset of oscillating models that 
has not been previously explored. In these models, the 
field is initially slowly rolling down the potential, with 
increasing w$; it then transitions into rapid oscillations 
at late times. An interesting feature in this regime is the 
non-smooth structure of the confidence contours. This is 
a real effect. 

To illustrate the reason for the non-smooth confidence 
contours, we have sampled four representative models in 
this regime, denoted by the four diamonds in Fig. [31 
which correspond to n — 0.16 and <f>i =0.10, 0.12, 0.13, 
and 0.15. (These are the four models previously displayed 
in Fig. [5]). This illustrates the complex structure of the 
allowed region as 4n = 0.10 and 0.13 are ruled out at 
95% confidence, while <pi = 0.12 and 0.15 are not. Next, 
we display the distance modulus for these four models, 
along with the supernova data, in Fig. [51 Although the 
individual models cannot be distinguished here, it is clear 
that all of these models lie near the edge of the allowed 
region because they tend to produce a distance modulus 
that is lower than the observations. In Fig. [6l we show 
a blow-up of Fig. [5] The smallest values of the distance 
modulus are produced for = 0.10 and 0.13; these are 
the two models that are ruled out at 95% confidence. The 
important point here is that the distance modulus is not 
a monotonic function of <f>i . The reason for this is illus- 
trated in Fig. where we give the redshift dependence 
of the density of the scalar field (which determines H(z), 
and therefore, ^{z)). Comparing this figure with Fig. 
[2 we see that H(z) depends crucially on the oscillation 
phase at the present day. The two models for which the 
field is at a maximum of the potential at present (thus 
giving = —1 at present) are 4>i = 0.10 and (f>i — 0.13. 
These produce a larger density at moderate redshift and 
therefore a smaller integrated value of 1/H(z). Thus, the 
complex behavior of our contours in Figs. [3] and 2] is a 
phase effect: the phase of (j> at present is an oscillating 
function of 4>i , and this phase determines whether or not 
a given model is excluded. For early oscillation models 
the phase becomes irrelevant, since these models oscillate 
many times before the present, and the equation of state 
parameter takes on its period-averaged value. 



IV. CONCLUSIONS 

Our results indicate that oscillating quintessence mod- 
els, even for very simple potentials, can display more 
complex behavior than has previously been considered. 
Previous discussions have centered on the early oscilla- 
tion case, in which the oscillation frequency is always 
much larger than the expansion rate when the scalar field 
dominates the expansion. Our results show that the late 
oscillation case, in which the scalar field transitions from 
a slowly rolling regime at early times into oscillations 
at late times, can also be consistent with current ob- 
servations, although the behavior of these latter models 
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FIG. 5: Distance modulus for the indicated values of <j>i along 
with the supernova values used in our analysis, for V(<f>) = 
*#| n , with n = 0.16. 
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FIG. 6: A blow-up of Fig. [5] illustrating the difference in 
behavior of the four indicated models. 



depends crucially on the oscillation phase of the scalar 
field at present. Johnson and Kamionkowksi [25[ have 
argued recently that rapidly-oscillating scalar fields with 
a time-averaged value of less than are unstable to 
the growth of inhomogeneities (see also the argument in 
Ref. [IB])- In this case, only the late oscillation models 
would remain as acceptable oscillating scalar field models 
for dark energy. 



How general are our results? The conclusion that v/H 
increases with time depends only on the assumption that 
the scalar field energy density decays more slowly than 
the background energy density, so it should be generally 
applicable to any plausible model for dark energy. On 
the other hand, our comparison of models with observa- 
tions was restricted to the case of power-law potentials. 
However, we would expect qualitatively similar results for 
other bound potentials with a single minimum at V — 0. 
If the minimum of the potential were at Vo > 0, the 
analysis is almost identical, except that the dark energy 
now has an additional component with constant density 
Vo- Such a model would have a values of w(z) closer to 
— 1 than the corresponding model with a minimum at 
V = 0, and thus, would tend to agree more closely with 
current observations. One type of potential that could 
display completely different behavior from those studied 
here is a V(4>) with multiple minima, such as the Mex- 
ican hat potential discussed in Ref. fl4| . In this case, 
the field evolution can be considerably more complex, as 
the field can initially oscillate over the entire range of the 
potential, but will eventually become trapped in a local 
minimum at late times. In this case the simple analysis 
presented here will not apply. 
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